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Abstract. The Pauli-Hamiltonian of a molecule with fixed nuclei in a strong 
constant magnetic field is asymptotic, in norm-resolvent sense, to an effective 
Hamiltonian which has the form of a multi-particle Schrodinger operator with 
interactions given by one-dimensional <5-potentials. We study this effective 
Hamiltonian in the case of the H2-molecule and establish existence of the 
ground state. We also show that the inter-nuclear equilibrium distance tends 
to as the field-strength tends to infinity. 

1. Introduction 

In the final decade before his untimely death, Pierre Duclos embarked upon 
a large-scale systematic study of atoms and molecules in strong homogeneous 
magnetic fields. This topic had already received considerable attention in the as- 
trophysics and mathematical physics community (see for example |AtMolj . [La] . 
[NeuSt] and [LSYj ) . The published record of Pierre's involvement started with the 
paper [BDlj (which can to a large measure be accredited to him), which established 
the norm-resolvent convergence, as the magnetic field strength tends to infinity, of 
the Pauli Hamiltonian of a non-relativistic atom in a constant magnetic field to 
an effective Hamiltonian describing a one-dimensional atom with the electrostatic 
interaction replaced by Dirac 8 point-interactions. Though there exists a long and 
respected tradition in theoretical physics and chemistry of using Hamiltonians with 
(5-potcntials as toy models (going back, at least, to the Kronig and Penney model 
of solid state physics) , this theorem gave an example of such (5-potentials naturally 
arising as part of a genuine physical problem. Around the same time, [BaSoYj 
proved by a different method that the ground-state energy of the true Hamilton- 
ian converges to that of the corresponding 5-Hamiltonian. The Hartree functional 
associated to the (5-Hamiltonian had already made its appearance in [LSY] . 

The approach initiated in [BD1] was further developed in [BD3] , which estab- 
lished a three-member hierarchy of effective multi-particle Hamiltonians on the 
line, each one of which is asymptotic, in norm-resolvent sense, to the full Pauli- 
Hamiltonian, but with increasingly better rates of convergence. The (5-Hamiltonian 
is the member of this family with the slowest convergence. This paper also an- 
alyzed the effect of imposing particle anti-symmetry which, on the level of the 
(5-Hamiltonian, translates into the latter acting on a certain direct sum of fcrmionic 
L 2 -spaces and L 2 -spaces with no symmetry-restriction at all; cf. section 2 below for 
the case of two electrons. An early announcement of some of the results of [BD3j . 
which took a long time in writing, was made in [BD2] . 

Despite its apparent simplicity, the spectral analysis of these (5-Hamiltonians is 
far from trivial, except in the case of a single electron. As a first step, Pierre, 
together with Santiago Perez-Oyarziin and two of the present authors, investigated 
in [B BDPj the molecular iJ^-ion in the Born-Oppenheimer approximation in strong 
magnetic fields, taking the corresponding (5-modcl as its starting point. The latter 
is explicitly solvable for a one electron molecule, and its solution can be used, in 
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combination with second order perturbation theory, to obtain the correct values of 
the equilibrium distance and of the binding energy for strong magnetic fields of the 
order of 10 9 to 10 14 . a.u. Building upon this paper, it was shown in |BBDPVT| (see 
also |BBDPV2j ) that as the field-strength B increases, the equilibrium distance 
between the nuclei in the Born-Oppenheimer approximation tends to 0, both for 
the 6-model and for the full molecular Pauli Hamiltonian, for the latter at a rate of 
(log_B)~ 3 / 2 . This naturally suggested the possibility of using strong magnetic fields 
to facilitate nuclear fusion by enhanced tunnelling through the Coulomb barrier, 
a point which was taken up by Ackcrmann and Hogrcve [AH3], who numerically 
computed accurate equilibrium separation distances, as well as the corresponding 
nuclear fusion rates, for a range of B's from and 10 4 . For field strengths B of the 
order of 10 4 they indeed found a drastic increase of the tunnelling cross sections. 
Unfortunately, such B's are still much too big to be generated in a laboratory on 
earth. 

In a parallel set of papers [CDRj and BBDR (see also |Bc| ) , as well as in hitherto 
unpublished work, Pierre, in cooperation with various co-authors and building upon 
an important earlier paper of Rosenthal [Ros] , developed a systematic method for 
analyzing the spectrum of these multi-particle <5-Hamiltonians, which he baptised 
the skeleton method, and which basically reduces the spectral analysis of such an 
Hamiltonian to that of a finite-dimensional system of (N — l)-dimensional integral 
operators, N being the number of particles (note that N — 1 is then the dimension of 
the support of any of the ^-potentials occurring in this 5-Hamiltonian) . This system, 
which is called the skeleton of the original Hamiltonian, is still too complicated to 
be solved explicitly even for an N as small as 2, but opens the way for a systematic 
numerical approach (as was done by Rosenthal for the Helium-type (5-atom with 
N = 2.). For arbitrary N the size of this system of integral operators will be N(N+ 
l)/2, though this can be reduced by particle symmetry and parity considerations. 

The present paper takes up the work of [BBDPVlj . [BBDPV2] for the H+- 
ion by studying the existence and equilibrium distance of the ^-molecule in the 
^-approximation. We show that the equilibrium distance, in this approximation, 
again tends to as B — > oo, at the same speed as that found for . The paper 
should be seen as a part of a larger and more ambitious project of Pierre, which 
was to prove this for arbitrary Pauli molecules in strong magnetic fields. In the 
context of enhanced fusion, one might hope, following | AH3| . that adding more 
electrons would lead to increased shrinking of the equilibrium distance through 
better shielding of the nuclear charges. Our results suggest that such an effect will 
not show up in the top order asymptotics of this equilibrium distance as function of 
the field strength B, though it might of course still make itself noticeable at finite 
but large B. A major problem for the analysis of the ^-molecule is that, unlike 
for H£ , the corresponding 6- model is no longer explicitly solvable, and we will 
arrive at our conclusions by a combination of analytical and numerical methods. 
In particular, a key intermediary result we need is that the electronic ground-state 
energy (that is, leaving out the internuclear repulsion term) of the molecular 6- 
Hamiltonian be an increasing function of the internuclear distance. This seems 
physically intuitive, but we have been unable to find an analytic proof. For non- 
relativistic single-electron molecules without magnetic field, monotonicity of the 
electronic energy in the internuclear distance was shown in [LSJ . [L] : we are not 
aware of any rigorous results for the multi-electron case, with or without magnetic 
field. By deriving explicit estimates for the equilibrium distance and minimal energy 
associated to an elementary variational upper bound of the molecular ground-state 
energy, and combining these with a trivial lower bound, we are able to show that it 
suffices to know the monotonicity of the (true) electronic ground-state for a certain 
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range of intcrnuclear distances. The latter was, amongst other things, verified 
numerically in [Be], using Pierre's skeleton method. 

As will be clear, this paper's conclusions are far from complete, and the paper 
is perhaps best seen as a work-in-progress report on a project Pierre was working 
on at the time of his death, one of the very many he was actively involved in. 
His absence, and the impossibility of the lively exchange of ideas we had grown 
accustomed to having with him, are sorely felt by the three surviving authors. 

2. Notation and Main Results 

As mentioned in the introduction, we want to study the equilibrium distance, 
in the Born-Oppcnheimer approximation, of the non-relativistic i^-molecule in 
a strong constant magnetic field, in the large field limit. We take the magnetic 
field directed along the z-axis, and the two nuclei, of charge Z, aligned in the 
direction of the fielcQ and located at ±\Re z , where e z = (0, 0, 1) is the unit vec- 
tor in the direction of the z-axis. The molecule is described by the familiar two- 
electron Pauli-Hamiltonian M B (which we won't write down explicitly) acting on 
anti-symmetric wave-functions which include spin-coordinates, and therefore are 
functions of (r,-,s,-), i = 1,2, where n are the spatial variables of the i-th elec- 
tron, and Si = ±1 are its the spin variables; the anti-symmetrisation is done with 
respect to all the variables. We fix the total (orbital) angular momentum in the 
field-direction to be M, which we can do since the corresponding component of the 
angular momentum operator commutes with the Pauli-Hamiltonian. In this situ- 
ation it can be shown by the methods of [BDlj . |BD3j that the full Hamiltonian 
ET 5 , after projection onto the lowest Landau band, is asymptotic, in norm-resolvent 
sense, to the following two-electron Schrodinger operator with 5-potentials on R 2 : 




Hs will act on a certain direct sum of copies of L 2 (R 2 ) and L\ s (R 2 ) which we 
specify below, where "a.s." stands for asymmetric wave- functions. The parameters 
a and e in the definition of Hs are related to the original parameters R, Z and B 
by 

(2) a := RLZ/2, e := Z/L, 

where L = L(B) = 2W(\/B /2), W being the principal branch of the Lambert- 
function, defined as that branch of the inverse of xe x which passes through and 
is positive for positive x; cf., |CGH JK] . Note that e = e(B) — > as B - >oo, since 
L(B) ~ logB. We also let 

s 

(3) h s := hs(a, Z) := H s - — , 

Za 

the electronic part of Hs (at fixed a). The operators Hs and hs are defined in 
form-sense, as closed quadratic forms on the first order Sobolev space H 1 (M 2 ). 
Their operator domain consists of those functions ip in H 1 (M. 2 ) whose restrictions 
to M 2 \ {zi — z 2 , Zi — ±a for i = 1,2} are in the second order Sobolev space H 2 , 
and whose gradients satisfy an appropriate jump-condition across the supports of 
the different ^-potentials: specifically d Zi 4>\ Zi =±a- — d Zi ip\ Zi= ± a + = 2tjj\ Zi= ± a . with 
a similar condition involving the normal derivative of tjj across {z\ = z 2 }, but with 
opposite sign: cf., the Appendix of }BD3j for details. 

^The methods of IBD1I . |BD3I on which this paper is based unfortunately do not apply if the 
molecule is not aligned with the field, since the angular momentum in the field direction is not 
conserved anymore. 
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We refer to [BD3| for the precise technical sense in which U B and Hg are as- 
ymptotic, but note that this will imply that if Hg has an eigenvalue E(L) < then 
M B will have an eigenvalue at distance O(L) of L 2 Z 2 E(L). 

The Hilbert space in which Hg and hg act is , by |BD3j . theorem 1.8, 
(4) L 2 (R 2 )* M '®Ll s XR 2 )* M2 , 

where Mi = {(mi, m 2 ) € N x N : mi + m 2 = M, mi < m 2 } and M.2 — {(m, m) : 
rn G N, 2m = M}, which is either empty or a singleton (depending on whether 
M is odd or even). Observe that unless M ^ 0, M.\ will always be non-empty, 
and, consequently, the infimum of the spectrum of Hg will be the infimum of its 
spectrum on L 2 (M. 2 ), which equals the infimum of the spectrum of its restriction 
to the bosonic subspace of symmetric wave- functions in L 2 (R 2 ), despite having 
originally started off with fermionic electrons. For the higher lying eigenvalues, no 
symmetry restrictions have to be taken into account, at least for the J-Hamiltonian 
- the situation will be different for the full Pauli-Hamiltonian. 

Since for M = 0, Hg acts on L 2 S (M?) while for M > 1, it acts on a Hilbert- 
space containing L 2 (R 2 ) as a direct summand, it is clear, for any fixed value of the 
inter-nuclear distance a, that the ground state of hg, if it exists, occurs in any of 
the sectors with M > 1. It then follows that for sufficiently large B, the ground 
state of (the electronic part of) the Pauli-Hamiltonian H s itself also exists and 
occurs in one of the sectors with M > 1. This was observed in numerical studies. 
The ground-state of the Hi molecule was computed numerically (by the variational 
method) in }DSC| . |DSDCj . and that of the two-electron He 2+ -ion more recently 
in |TGj . It was shown that for both these molecules the ground-state changes, 
with increasing B, from a spin singlet state with M = (for B = 0) to a, bound 
or unbound, spin triplet state (both electron spins parallel) still with M = to a 
strongly bound spin triplet state with M = 1 (using molecular term symbols, the 
transition is 1 S„ — > 3 S U — > 3 II U ). A similar phenomenon occurs for other small 
molecular ions such as H^ and HeH + , cf. |Tuj and its references. We should 
note that it is implicit in the definition of Hg that all the electron spins are taken 
anti-parallel to the magnetic field (since Hg effectively operates on the projection 
of the full Hilbert space with spin onto the lowest Landau state - cf. |BD3| ) so that 
the ground state of Hg will correspond to a triplet state of the Pauli-Hamiltonian. 
For further numerical studies of spectra of atoms and molecules in strong magnetic 
fields, sec for example the two conference proceedings [AtMolj and [NeuStj . 

Let -ff 2 (<5) denote the one-dimensional -ff 2 -molecule as defined by the <5-Hamiltonian, 
Hg. We will study this molecule for arbitrary nuclear charge Z, but will often single 
out the case of Z — 1, corresponding to ff 2 , as well as the case of Z = 2 which pro- 
vides an aymptotic description of the molecular i?e 2 + -ion in the large field limit. 

The two main results of our paper are: 
Theorem 2.1. (Existence of H 2(8)) 

(i) If Z > 1, then the electronic Hamiltonian hg possesses a ground state for all 
a > 0. 

(ii) The 77 2 ((5) -molecule exists (in Born Oppenheimer sense) for all Z > 1, as long 
as Z/L < 0.297. For Z = 2, this can be sharpened to 2/L < 0.458 or L^ 1 < 0.229. 

We recall that the molecule exists in Born-Oppcnheimer sense if inf a E{a) < 
lim inf a-j.oo E(a), where E(a) is the infimum of the spectrum of Hg = Hg(a), and 
if moreover there exists an a = a cq such that E(a cq ) = inf a E(a) and E(a cq ) is 
an eigenvalue of Hg(a eq ); note that inf a E{a) < liminfa-^x, E{a) and continuity of 
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E(a) already imply that the inhmum is attained. We will call any a eq in which 
E(a) assumes its global minimum an equilibrium distance of the molecule. 

Theorem 2.2. For fixed e, let a oq (e) be an equilibrium distance of the H2{5)- 
molecule. Then there exists a constant c = c(Z) such that 

(5) a cq (e) ~ c-v/e, as e — > 0. 

The constant is the same for any of the, potentially multiple, equilibrium dis- 
tances; it is natural to conjecture that a eq (e) is in fact unique, as is observed 
numerically: see for example figure |4] below. 

We conjecture that, as a corollary of these theorems, the actual molecule as 
modelled by the Pauli Hamiltonian can also be shown rigorously to exist for Z > 1 
and sufficiently large B, and that its equilibrium distance behaves as C(log\B)~ 3 / 2 
for some constant C (note that by ©, R ~ L~ 3//2 if a ~ C\fe). For this was 
shown to be the case in BBDPV1 . We recall that for B = and Z = 1, the 
existence of the -f^-molecule was rigorously established in |RFGSi . A discussion 
of the stability scenario for varying Z, still with vanishing magnetic field, can be 
found in [AHlj . |AH2j . where it is in particular shown that He\ + is meta-stable 
when B = 0. We also note that reduction of the size of the molecule with increasing 
magnetic field strength has been observed numerically: cf. for example |Tuj and its 
references. 



3. Existence of the molecule 

The ground-state energy of the one-electron Hamiltonian hg := — ^^-z+J2± 
a) on L 2 (R), which is defining the electronic energy of the (<5)-ion, was deter- 
mined in [BBDPj as being — ^ao(a) 2 , where 

W(2ae- 2a ) 

(6) a (a):=l + ^^ L , 

with normalized eigenfunction <po given by 

A ie - ao ^, \z\>a 



(7) M*) = 

where 



A 2 cosh(aoz), \z\ < a, 



(8) Ai = , ^( 1 + e2a "°) A2 - 



1 + 2e 2a "° + 2aa ) + 2e 2aa ° + 2aa ) 

In fact, for a's larger than some ^-dependent threshold, the spectrum of hg 1 ^ con- 
sists of two eigenvalues ao(a), ct\ (a) less than and a continuous spectrum equal to 
[0, oo); for a's below this threshold the excited state is absorbed in the continuous 
spectrum, cf. |AM| and jHo| . The following properties of ao(a) will be repeatedly 
used below: ao(a) is a decreasing function of a, Qo(0) = 2 and ao(a) — > 1 as a — > oo. 
Finally, 

(9) a (a) - 2-4a + 0(a 2 ), a -> 0. 



Recall that the two-electron Hamiltonian hg = hs(a,Z) defined by © corre- 
sponds to the electronic part of the energy of our ^(o^-molecule. Let e(a) = e(a, Z) 
be the infimum of the spectrum of hg(a, Z), and 

(10) E{a) := E(a, Z, e) := e(a, Z) + ±, 
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the ground-state energy of the molecule at fixed internuclear distancfl To prove 
that e(a) is an eigenvalue we have to show, according to the HVZ-theorem, that 
e(a) is strictly less than the inf of the spectrum of that is: 

(11) e(a) < -ia (a) 2 . 

We can derive a simple variational upper bound for e{a) by using ipo ® ipo as test- 
function. Explicit evaluation of 

(5(zi~ Z2),(pl(z 1 )ipl(z 2 )) = / <p (z) A dz, 

JR. 

then leads to the upper bound 

(12) e X3B (a,Z):=-a (a) 2 + Z- 1 f(a) 
where 

, . . . 8cosh 4 (aao) + sinh(4aao) + 8 sinh(2aao) + 12aao 

(13) /(a) = ao 4( e »-o + 2aa + l) 2 ' 

which can also be expressed as 

e iaa ° + Ae 2aaa + 4sinh(2aa ) + 12aa + 3 

(14) /(a)=a ° 4( e ^o +2aao + 1)2 ■ 

It is elementary to show that /(a) < ^ao(a), the inequality being sharp for 
a = 0. The simple upper bound (|12p suffices to show that the ground state of hs(a) 
exists for any a > 0: 

Proof of theoremWIYi). It suffices to show that e UB (a, Z) < -\a Q (a) 2 for Z > 1. 
But this follows from e UB (a, Z) < -a (a) 2 + a Q {a)/2Z < -a {a) 2 + a (o)/2 < 
— ao(a) 2 /2, since ao(a) > 1 for a > 0. □ 

To show existence of the molecule, we examine when min a E VB (a) < lim inf _k>o E(a) 
where 

(15) E VB (a) := £ UB (a, Z, e) := e UB (a) + ^, 

the variational upper bound for the molecule's energy. We note the trivial lower 
bound: 

(16) E(a, Z, e) > E(a, oo, e) = -a (a) 2 + ^- =: E m (a 7 e), 

2a 

where "NI" stands for non-interacting electrons. 

Clearly, liminf a _ ) . 00 E(a) > lima-joo i? NI (a) = lim a _i. 0O — «o(a) 2 = — lj so to show 
existence of the molecule it suffices to show that there exists an a > such that 
E VB (a, Z,e) < —1, or, equivalently, e < j(a, Z), where 

(17) j(a, Z) := 2a (a (a) 2 - 1 - Z^fta)) . 

To derive a result which is valid for all Z > 1, we note that Z > 1 implies that 
j(a, Z) > j(a, 1). The molecule will therefore exist for any value of the parameter 
e = Z/L which is less that max a j(a, 1), for any Z > 1. For the case of Z = 2, we 
have to have that e < max a j(a,2). Part (ii) of theorem 12.11 is then an immediate 
consequence of the following lemma: 

Lemma 3.1. j{a, 1) has a global maximum of 0.297 on a > 0, which is attained for 
a = 0.254. The maximum of j(a, 2) is equal to 0.458, and is attained for a = 0.337. 



2 as a notational convention, we denote electronic Hamiltonians and energies by lower case 
letters h$, e(a), e UB (a), etc. and the corresponding molecular entities (obtained by adding e/2a) 
by upper case letters Hg, E(a), B UB (a), etc. 
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This lemma has been checked numerically, using Maple: see figure Q] for the 
graphs of j(a, 1) and j(a, 2). 




FIGURE 1. Graphs of j(a, 1) ( — ) and j( a ,2) ( ) 



Remarks 3.2. (i) Our reliance, here and below, on numerical analysis may be felt 
to be somewhat unsatisfactory, and it may be possible, with sufficient effort, to 
give a rigorous analytic proof of lemma I3~T1 since j(a, Z) has an explicit analytic 
expression (and similarly for g[a,Z) below). This lemma of course concerns only 
one particular variational bound (using a simple test-function) which is unlikely to 
be optimal, although numerically it does not perform badly (cf. figure 2] below). It 
would also be interesting to find better variational bounds, if possible again with 
closed analytic expressions. 

(ii) The Taylor expansion of j(a, Z) in a = starts of as j(a, Z)/2 = (3 — Z~ 1 )a + 
0(a 2 ), so that max a j(a, Z) will be strictly positive if Z > 1/3. Graphical analysis 
shows that j(a,Z) < when Z = 1/3 and therefore for all Z < 1/3 (note that 
linia-i.oo j (a, Z) = — oo). These remarks imply that the present method would prove 
existence of the molecule for all Z > 1/3, for sufficiently large L(B), provided we 
can show that hg(a, Z) possesses a ground state for such Z, at least for those a for 
which j(a, Z) > 0. 

To determine min a i? UB (a, Z), we have to solve <9 e UB = e/2a 2 , or 
(18) g(a,Z)=E, 

where g(a,Z) := 2a 2 (-2a (a)a' (a) + Z' 1 f'{a)), the prime indicating differenti- 
ation. Graphical analysis of g(a, Z) using Maple or Mathematica shows that (JTSJ) 
has two solutions as long as e < max a g{a, Z) - sec figure [5] below. 
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Figure 2. Graphs of g(a, 1) ( — ) and g(a, 2) ( ) 

It follows (again only graphically, for the time being) that E VB (a, Z 7 e) has two 
critical points, of which the smaller one, for sufficiently small e, turns out to be a 
local minimum and the larger a local maximum^. The local minimum is of course 
a global one under the hypotheses of theorem 12. lf ii). 

Let the minimum of £' UB (a, Z, e) be attained in a = a^ B (e); a u q B (e) would be 
the equilibrium position of the molecule if its energy were in fact equal to our upper 
bound. It then follows (from the graph of g(a, Z)) that et^ B (e) — > as e — > 0. This 
information turns out to be sufficient to determine its small-e asymptotics, using a 
simple argument which we will also use in the next section for the true equilibrium 
position of iJ2(<5)-molccule: 

Lemma 3.3. As e -> 0, 

(19) a UB ( £ )^ C Vi, 
with constant c equal to 

1 _ 1 / Z 

( ' °~ y/2(e™)>ffi ~ 2 V 8Z- 1" 

Proof. We fix Z and write e UB (a, Z) = e UB (a). Then a™(e) is the smallest of the 
two solutions, a = a_(e), of (e UB )'(a) = (2a 2 ) _1 e, so that 

a - {£) = \[~2(e VB na-(e)Y 

^as follows by further numerical analysis: for example, in the stationary points of B UB (a), we 
have that (E UB )"(a) = -2a' Q {o) 2 -2a (a)a'^(a) + Z- 1 f" (a) + a- 3 g(a, Z), and the latter function 
is positive for a in a neighborhood of (e.g. for a < 2 in case Z = 2) 
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Now, as we observed above, a_(e) — > as e — > 0, and therefore (e UB )'(a_ (e)) — > 
(e UB )'(0), from which the first asymptotic equality of (fH?|) follows. For the second 
equality, we use that e UB (a) = -(4 - Z" 1 ) + (16 - 2Z- 1 )a + 0{a 2 ). □ 

Note that the argument we gave here is very general, and we will seek to apply 
it in the next section to the true equilibrium position, a eq (e) of H2(S). 

The equilibrium position a^ B (e) we just studied has of course no real physical 
significance. We will use it, though, in the next section, to obtain a weak a-priori 
estimate for a oq (e). To that effect, we note that if we combine (fH)]) with the Taylor 
expansion of E VB (a,e) we find: 

Corollary 3.4. There exists a constant C > such that 

(21) ^fn(e) :=min£; UB (a,e)~-4+Z- 1 +C x /i, e -> 0. 

a 

In particular, E^f n (e) < — 1 for Z > 1 and £ sufficiently small . 
One can perform a similar analysis for the lower bound, 

E m (a,e) = ~a (a) 2 + ^- , 

and show that the curve -E NI (a,e) has, for fixed and small enough e, a local mini- 
mum and a local maximum, the local minimum again being absolute for sufficiently 
small e; the local maximum lies above lim^oo E m (a,e) = —1. If a^e) is the lo- 
cation of the minimum, then we have again that a^q(£) — Cyfi with a constant c 
which is now equal to c = (2(e NI )'(0)) -1 / 2 , and 

(22) min £ NI (a, e) - -4 + Cy/s, 

a 

for a suitable constant C. We finally observe, from the geometry of the graph of 
E m (-,e), that E NI (a) will be strictly increasing on any interval [a^J(e),A] such 
that E m {A,e) < -1 = Iim^oo E m (a). 

4. Asymptotics of a eq (e) 

Recall that e(a) := e(a, Z) is the ground-state energy of h$, and that E{a 1 e) = 
e{a) +e/2a. As we have seen, E(a, e) possesses an absolute minimum if e = Z/L is 
sufficiently small. It is not a priori known whether this minimum is unique, though 
we would expect it to be. Below, we let a oq (e) be any value of a at which E(a,e) 
attains its absolute minimum in a on [0,oo). 

We will use the following lemma on differentiability of e(a), whose proof we 
postpone till the the end of this section, in order not to interrupt the flow of the 
argument. 

Lemma 4.1. The ground-state energy e(a) is continuously differ entiable on [0, oo). 
In particular, its right- derivative e'(0 + ) in exists. Moreover, e'(0 + ) > for 
Z > 1/4. 

The next lemma takes up the idea of lemma [3~B1 but now for e(a) directly. 

Lemma 4.2. Let a = a(e) be a solution of e'(a) = e/2a 2 for which there exists a 
constant A > such that 

(1) a(e) S [0, A] for s sufficiently small; 

(2) min oe[ o,A] e'(a) > 0. 
Then as e — > 0, 

(23) a cq (e) * ^^y- 
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Proof. Since 2a(e) 2 = e/e'(a(e)) < e/minmyt] e'(a) — > as e — > 0, it follows that 
e'(o(e)) -> e'(0+), and hence that 2a(e) 2 /e e'(0+)" 1 . □ 

We next use E m (a,e) and E^f n (e) (cf. lemma I3T41 to effectively bound a eq (e) 
(compare [BBDPVTj . proof of theorem 2). We assume e sufficiently small for 
E VB (e) < -1 to hold: cf. corollary 1X1 

Lemma 4.3. Suppose that e is sufficiently small so that E^f n (e) < — 1. Let A = 
A+(e) be the largest of the two roots of E m (A,e) = E^f n (e). Then a cq (e) < A + (e). 

Proof. We first claim that E m (a cq (e), e) < — 1. For otherwise, E(a cq (e),e) > 
E NI (a cq (e),e) > -1 which is in contradiction with E UB (e) < -1 (since E^(e) 
dominates the ground-state energy of Hs). Suppose now that a eq (e) is strictly 
larger than A + (e) which, as the larger root, is bigger than a^(e). It then follows 
from the properties of a — > E m (a, e) that E m (a, e) will be strictly increasing on the 
interval [a^ q (e),a eq (e)], and hence E(a eq (e),e) > E m (a eq (e), e) > E m (A + (e), e) = 
-^mki( e ) — E(a cq (e) 7 e), which is a contradiction. □ 

One might hope that the previous lemma, combined with lemma 13.31 would 
already imply that a eq (e) — > 0, in which case lemma B~2l together with e'(0+) > 
would already imply theorem 12.21 However, this is unfortunately not the case: if 
A + (e) -> 0, then since A + (e) > fl™(e) - yfi, 

E m {A + {e), e) < a (A + (e)) 2 + c^e -> a (0) 2 = -4, 

but on the other hand E m (A + (e),e) = E^f n (e) ->• -3 + Z" 1 , by (J2TJ, which is a 
contradiction. However, it is easy to show that a eq (e) = 0(1) as e — > 0, since if 
A + {e v ) — > oo on some sequence e v — > 0, then 

E m (A + (s„) 7 s„) = -av{A + {e v )) 2 + h- -1, 

which contradicts 

E m (A + (e„), e v ) = E^f n (e v ) ~ -4 + Z" 1 + -> -4 + Z"\ 

as long as Z > 1/3. // we now would know that e'(a) > /or all a > 0, this 0(1)- 
estimate for a cq (e) in combination with lemma |4~21 would prove theorem 12.21 We 
were however not able to prove monotonicity of e(a). What we can do is give an 
effective upper bound for a eq (e), which means we only have to check monotonicity 
numerically on some known finite interval. 

Lemma 4.4. Let Z > 1 and let e be such that E^f n (e) < —2: (this is true for all 
sufficiently small e by Corollary \3.4\ ). Then A + (e) < o.q (v^) = 0.3116. Conse- 
quently, a eq (e) < 0.3116. 

Proof. Elementary, if we draw the graphs of E NI (a,e) = — a^(a) 2 + e/2a, of 
— ao( a ) 2 j which is an increasing function from -4 in to -1 at infinity, and of the 
constant functions E™ n (e) and —2. □ 

We can therefore take 1 (v / 2) in lemma H?2] (smaller A's are also possible, e.g. 
for Z = 2, depending on how large we are willing to let Z be or how small e). To 
finish the proof of thcorcm l2.2l we verify numerically that e(a) is strictly increasing 
on [0,a^\V2)}. This can be done using Pierre Duclos' skeleton method to compute 
e(a), and has been carried out by one of us in |Bej . The plot of e(a) reveals its 
monotonicity over the desired interval, thereby completing the proof of theorem l2.2l 
see figure[31 which for comparison also shows the one-electron energy — iao(a) 2 , as 
well as e UB (a) and e NI (a). We note in particular, here and also in figure 4 below, 
that our variational upper bound, though relatively naive, gives a quite reasonable 
approximation to the actual ground-state energy. 
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Figure 3. Electronic energy curves; from top to bottom: 
-iao(a) 2 (- - ); e UB (a) (---); e(o) (- -); e NI (a) (- - - ) 



The skeleton method allows us to numerically compute the equilibrium distance 
and energy of the H 2 (<5)-molecule for a given Z and L. We discuss by way of 
example the case of a Hydrogen molecule (Z = 1) in a magnetic field with field- 
strength corresponding to L = 10, or B = L 2 e L ~ 2.2 10 6 . We recall (cf. for 

3 3 

example [BBDP]) that the magnetic field is measured in units of Bq = m || = 
2.35 x 10 9 G = 2.35 x 10 5 T, where G stands for Gauss and T for Tcsla. Moreover, 
the energy is measure in units of hojQ = h^ 3 - = 27.2eV = 1 Hartrcc, where u>o is the 

cyclotron frequency, and distance in units of the Bohr radius ao = 2 = 0.53A 
Our L = 10 therefore corresponds to a magnetic field of 5.17 10 n T, which is of 
course not realizable on earth, but may be realistic for a neutron star. Figure U] 
displays the graph E(a) (solid line) as well as, for the sake of comparison, that of 
E VB (a) (medium dash line), E m (a) (small dash line) and — ^ceo{a) 2 + e/2a (large 
dash line). We see that a eq w 0.1 and E eq as —1.75, which corresponds to an 
equilibrium energy of hioo X —1.75 = —A7.6eV = — 1.75Hartree and an equilibrium 
distance of R eq = ao x as 10 _2 A. The equilibrium distance is much smaller 
than ao but still significantly bigger than the distance of as 10~ 5 A over which the 
nuclear interaction between two protons makes itself felt. Nevertheless, following 
[AH3] . it may be small enough to significantly enhance the probability for protons 
to pass through the electronic barrier and be trapped in the nuclear well. It would 
be interesting to compute the tunnelling cross-section in the Gamov model for the 
i?2 (<5)-model, as simplified model for the actual H2 molecule. 

We finish with the promised proof of lemma 14.11 

Proof of lemma \4-l\ We will carry out the proof for the iV-electron Hamiltonian 
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where we will suppose, to fix ideas, that a > (though this is not strictly speaking 
necessary). Differentiability of e(a) follows from norm-differentiability of the resol- 
vent, which can be established using the symmetrized resolvent equation - we skip 
the details. 

We next study e'(a). Let ip a be the normalized ground-state eigenfunction of 
hs(a) on L 2 (M. N ), which is unique. A formal application of the Feynman-Hcllman 
theorem would lead to 

e '( fl ) = - E E (*'(* ± °)^. vv) = - E E ± a )> i^i 2 )- 
i ± i ± 

However, we have to be careful here: since d Zi ip a a priori has a jump in Zi = ±a, 
\ip a \ 2 is not an admissible test- function for S'(zi ± a). 

We will first establish Feynman-Hellman on the level of quadratic forms. If we 
write ipa{z) = tp(z,a) (to more clearly bring out the parameter dependence in the 
notations), the eigen-equation hs(a)tp a = e(a)ip a translates into 

J(VV>G,a),V^)-X]E / / ^(•,o)? = e(o)(^(-,o),VJ), 

for all ip in iJ 1 ^), the form-domain of h s (a) We now carefully differentiate with 
respect to a, taking into account that ifj{-,a) is only left- and right- differcntiablc 
with respect to Zi in Zi = ±o. Furthermore, the derivative of the arbitrary H 1 - 
function ip only exists in L 2 -sense. However, anticipating that we will take <p = 
ip a = ip(-, a) below, we assume that the relevant left- and right- derivatives of <p also 
exist. Using the norm-differentiability of ip(-,a) with respect to a (as an L 2 -valued 
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function of a) and taking right-derivatives with respect to a, we arrive at: 

-(Vd a ip(-,a),\7ip) - y^y~] [ d a ip(zi,...±a,...,ZN,a)lp(zi,...,±a,...ZN) 

1 i ± JR"- 1 





f d Zi i/;(zi,.. 

R»-i 


. a + , . . . ,zjv,a)^(2i, . . . ,a, . . . 






1 i>(zi, . . . a, 

RJV-1 


•■•,ZAr,a)9 Zi ^(zi,...,a + ,... 


,zn) 


+ ?j 


f d Zi i/}(zi,.. 


. - a - , . . . ,z w ,a)^(zi, . . . ,a, 


• • • , 2jv) 


+ ?j 


f ^(zi,...a, 

JJJV-l 


. ..,Zjv,a)e? Zl p(zi,...,-a _ ,. 


..,2jy) 



= e(a)(d a ip(-,a), <p) + e'(a)(ip(-, a),<p). 

(Note that the right-derivative (with respect to a) of a term such as ip(zi, . . . , —a, . . . zn) 
is —d Zi ip{z\, . . . , — a - , . . . , Zjv).) As announced, we now take <p = ^(-j a). Writing 
again ip a for ip(-,a) and using real-valuedness of -0a (being the ground state func- 
tion), we then arrive at the identity 



{hsiajdaifra^a) - V" / 2^ a {z\ , ■ . ■ , a, . . . , z N ) d Zi ip a (zi , . . . a + , . . . , z N ) 

i it"- 1 

+ 2tp a (zi,...,-a,...,Z N )d Zi ip a (zi,...-a~,...,Z N ) 

= e{a){d a ^j a ,ipa)+e!{a){^ a ,ijj a ). 

Since (hs(a)d a ->pa,ipa) = {d a ^a,h s (a)-)p a ) = e(a)(d a ipa,il>a), we have proved: 

Lemma 4.5. (Feynman-Hellman for h$(a)). Letip a be the normalized ground- state 
ofh(a). Then 



(25) 



End of the proof of lemma \^T\ If we now let a —> in (|25p. and use the boundary 
conditions at a = for membership of the domain of h$(§), we obtain that 

e'(0+) = 2 V f (d z ^oL=o- ~ a Zi | 2!=0+ )^o 

i J*i=0 

= / O 2 >O. 

Since, finally, 

2 E / = i||V^o|| 2 +E( (5 ^-^)V'o,0o)-(^(O)Vo,V'o) 

i 1 i<j 

> -(h s (0)ih,1>o) = -e(0), 

we find, specializing to TV = 2 again, that e'(0+) > -2e(0) > -2e UB (0) = 2a (0) 2 - 
2Z- 1 f{0) = 2a Q (0) 2 - Z- l a Q (Q) > as long as Z > (2a {0))- 1 = 1/4. □ 
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5. Conclusions 

Wc have approximated the ^-molecule in a constant magnetic field, as described 
by the non-relativistic Pauli Hamiltonian with fixed nuclei, by a two-electron model- 
Harniltonian of a one-dimensional molecule with electron-electron and electron- 
nuclei interaction given by (5-potentials, and interaction between the two nuclei 
given by the usual Coulomb potential. It can be shown, using the methods of 
[BD3] . that this approximation is exact in the large field limit. We have shown, for 
this approximation, that the ground state of the molecule exists, for any nuclear 
charge Z > 1 and for sufficiently large magnetic fields, and that the (re-scaled) inter- 
nuclear equilibrium distance of the molecule tends to with increasing field-strength 
B, at a rate of (log B)^ 1 / 2 . This generalises earlier results for in [BBDPj . 
[BBDPVl] . [BBDPV2 . Of the numerous questions which remain, the foremost 
one is to extend these results to the full Pauli-Hamiltonian. This was possible for 
the single-electron 77^"-ion, but the argument for the equilibrium distance given in 
[BBDPVT] used the exact solution of the 5-model. In the case of the two-electron 
Hi such exact solutions are not available anymore, despite the apparent simplicity 
of the (5-potentials, and we have had in part to rely on numerical computations to 
prove our results, notably to establish monotonicity of the electronic ground-state 
energy as function of the intcr-nuclcar distance. Since this would already imply 
that the inter-nuclear distance tends to at the proper rate (cf. the remarks just 
before lemma PO| . it would be very interesting to find an analytical proof of this 
monotonicity, or at least on a sufficiently large interval (such as the one specified in 
lemma EOT| . Another interesting question which remains open is that of uniqueness 
of the equilibrium position, for the ^-approximation as well as for the full Pauli- 
molecule. 

We also relied on numerical computations in the study of our variational upper 
bound and lower bounds. Here it may be possible to give analytic proofs, since 
these upper and lower bounds have analytically closed expressions involving a par- 
ticular and well-understood special function, the Lambert VF-function. It would 
furthermore be interesting to find sharper variational upper bounds, e.g. by using 
test functions which incorporate electron-electron correlations. 

Other perspectives are to go beyond the Born-Oppcnheimcr approximation and 
analyse the effect of nuclear vibrations on the stability of the nuclei, and to com- 
pute the probability of tunnelling through the Coulomb barrier. Finally, it would 
be interesting to see to what extend the methods and results of this paper can 
be extended to diatomic molecules with an arbitrary number of electrons. Some 
(unpublished) work in this direction was already started by Pierre Duclos. 
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